Multipartite transmission of quantum solitons 
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We analyze the preparation (launching) of a ID-propagating multipartite quantum soliton, its scat- 
tering and disintegration by an external potential. The regimes of suppressed disintegration and 
atom-number-dependent transmission are identified. 
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Thus far, most of the work related to solitons has 
dealt with their nonlinear wave properties [1], along with 
quantum noise (squeezing) corrections to their mean- 
amplitude ( "mean- field" dynamics [2,3]. However, the 
fascinating effects of second-quantized nonlinear dynam- 
ics that are known to date [2-6] suggest that a much 
richer variety of novel non-mean-field properties, partic- 
ularly quantum soliton features, is in store in nonlinear 
many-body systems. 

The quest for quantized solitons may be regarded as 
part of a broader effort to understand multipartite en- 
tanglement in complex quantum systems that have so 
far been analyzed at the mean-field level, such as BECs 
in traps [7] or optical lattices [8-11], whose dynamics is 
highlighted by Mott insulator transitions [9] and matter 
gap solitons [10,11]. 

Perhaps the most powerful tool of the non-mean-field 
many-body theory is the Bethe ansatz, ascending to 
Bethe's early work [12] and providing an exact quantum- 
mechanical expression for a one-dimensional (ID) quan- 
tum soliton, consisting of pairwise interacting N particles 
[13]. This solution has been extended to a two-band pe- 
riodic structure with Kerr nonlinearity, where band-gap 
quantum solitons have been predicted [14]. 

We analyze the preparation (launching) of a ID- 
propagating multipartite quantum soliton, its scatter- 
ing and disintegration by an external potential. The 
regimes of suppressed disintegration and atom-number- 
dependent transmission are identified. 

Let us consider a ID system of N identical bosonic 
particles described by the Hamiltonian 
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d 2 



2m e ff dx 2 
&(x)¥(x)i>(x)if(x)\ , 
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where $(x) is the bosonic-field annihilation operator, 
m e ff is the effective mass of the particle (which, in the 
presence of a lattice potential, may substantially differ 
from the free particle mass), g is the coupling constant 
of their ID contact interaction, and U(x) is a localized 



external potential. The effect of a possibly present peri- 
odic lattice potential has been already taken into account 
by introducing the effective mass. 

Exact eigenstates of the Hamiltonian (1) are given in 
the case of vanishing external potential, U(x) = 0, by the 
Bethe ansatz [12-14]: 
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where the function fa (x\, ... ,Xn) is a stationary solution 
of the linear TV-particle Schrodinger equation 
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An TV-particle bound state [15] requires effective attrac- 
tion, m e f f g I (2h 2 ) = — k < 0, i.e., either the scattering 
length or the effective mass (determined by the band cur- 
vature in a lattice) must be negative. This bound state 
has the wavefunction 



fx\xi, ...,x N ) 
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(TV-1)!(2k)( w - 1 )/ 2 

IpbN = Tttt ex P 



-K^2\Xj -Xi\ J , 
3<l 



(5) 

where X = {x\ + ... + x N )/N and NUK are the 
co-ordinate and momentum of the system's center- 
of-mass (CM), respectively. The coefficient on the 
right-hand side of Eq. (5) ensures the normalization 
J dx\ ...J dxjv|/jv^| 2 = 1j where the integrals are taken 
over the quantization length £ — > oo. 

The eigenenergy corresponding to the state (5) is 



e b (N) = —N(N - l)h z K 2 /(6m eff ). 
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The N 3 — N w iV 3 dependence of Sb(N) may allow the 
discrimination of a Fock state \N) if the difference be- 
tween Sb(N) and Eb(N — 1) is spectroscopically resolv- 
able. 

In the TV 3> 1 limit we may take Xj— X to be statitically 
independent for different j's, their distribution given by 
the square of the hyperbolic-secant solution of the clas- 
sical nonlinear Schrodinger equation [1] with the mean 
square value (xj — X) 2 — j^(Nk)~ 2 . Practically, the 
classical limit for the interatomic correlations is attained 
for N > 10 [3]. The distinct feature of the bound state 
(4, 5) is that the interparticle correlations are immune to 
dispersion. 

Let us now consider a quantum soliton incident on an 
external potential U(x) and analyze the quantum soliton 
disintegration probabiliy as compared to its transmission 
or reflection probability. To this end, we represent the 
potential acting on each of the N particles in the quan- 
tum soliton as a Taylor series 

N N 

]T U( Xj ) = NU(X) + -U"(X) J2(*j X f + ... . (7) 
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If the r.m.s. deviation of the jth particle's position from 
the CM is small compared to the characteristic length 
scale or 1 , over which U(xj) varies, then only two terms 
in the Taylor expansion (7) need be retained. The leading 
term, NU(X), then accounts for the potential felt by the 
bound multiparticlc complex that moves as a particle- like 
entity of mass Nm e ff, whereas the second (quadratic) 
term accounts for disintegration of that entity, due to 
the slight differences in the potentials felt by individual 
particles: an analog of a tidal force. The truncation of 
the Taylor expansion (7) is validated for 



N > ol/k. 



(8) 



As long as the disintegration probability is small, the 
bound state can be described by the wavefunction 



f N {xi, ...,X N ) = -j®NK(X)lpbN, 



(9) 



where -06 n is given by Eq. (5) and the CM motional 
state is the solution of the Schrodinger equation 
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We now choose the form of the external potential to 



be 



U(X) = -U cos\C 2 {aX), 



(11) 



because of its convenient properties and proximity to a 
Gaussian (see below). The CM wave function, which is 



expressed via the hypergeometric function, reduces to the 
sum of the incident and reflected waves at X — ► — oo and 
to the transmitted wave at X — > +oo. The corresponding 
transmission coefficient is found to be [16] 



sinh 2 (7riVisr/a) 
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A simple numerical check shows that similar transmis- 
sion dependence on the particle number of the slowly 
moving quantum soliton is attained also for a Gaussian 
potential, U(X) = -V exp(-a 2 X 2 ), if V w 1.56 U . 

We have calculated the bound-state disintegration 
probability V upon passsage through the potential us- 
ing the Fermi golden rule. This calculation admits three 
distinct regimes of the bound-state disintegration: 

(a) Disintegration without the change of the CM mo- 
tion. If \s f - s b \ ~ \de b /{dN)\ < Nh 2 K 2 /(2m e// ), 
then perturbations of the relative interparticle motion 
by the "tidal" forces do not affect the CM motion. In 
this case we arrive at the following estimation: V < 
(tt 6 / 180) N [a /{Nk)] 4 . This destruction probability is 
typically very small, provided N 3> (ch/k) 4 / 3 . Thus 
TV ~ (q;/k) 4 / 3 is the threshold for the disintegration to 
become appreciable. 

(b) Disintegration with CM acceleration. For 
\de b /(dN)\ > Nh 2 K 2 /{2\m eff \) and m eff < 0, a 
negative effective mass characterizes band-gap solitons 
near the bottom of the gap, the removal of an atom 
from the bound state by the "tidal" forces and placing 
it into the continuum is accompanied by energy release 
that accelerates the CM motion. If N 3> (a/n) 2 , the 
change of the CM momentum is sufficient to make the 
bound-state destruction probability exponentially small: 
V oc exp(— n^/Nn/a). Here N ~ (a/n) 2 is the disinte- 
gration threshold. 

(c) Disintegration with trapping by the potential. For 
\ds b /(dN)\ > Nh 2 K 2 /(2m eff ) and m eff > 0, the de- 
struction of the iV-particle complex is energy-consuming. 
As a result, the remaining (undestroycd) "core" with a 
reduced number of bound particles, N' < N. can be 
trapped in one of the discrete (bound) states of the po- 
tential N'U(X). However, if N > j(a/n) 2 , the released 
energy of internal motion becomes larger than the poten- 
tial well depth, NUo- Therefore TV <~ \{ajn) 2 is the 
disintegration threshold in this case. For N above this 
threshold, energy conservation prohibits in this case any 
change in both the internal and CM motion of the N- 
particle quantum soliton, so that V = 0. 

We now proceed to discuss an experimental scenario 
that may allow the observation of the effects described 
above. We consider ultracold alkaline atoms with the 
mass m in a wavegude, namely, a trap with strong trans- 
verse (radial) and weak longitudinal confinement. The 
atomic s-wave scattering length a can be enhanced by 
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means of a Feshbach resonance [17], up to \a\ s=s w r , A.J. Kerman, C. Chin, and S. Chu, ibid., 82, 1406 (1999); 

w r being the characteristic radial confinement. The ef- J- Stenger et al, ibid., 82, 2422 (1999). 

fective coupling constant of ID contact interactions is 
g = 2h 2 a/w^m 1 allowing one to achieve k ~ w" 1 . For 
laser-induced potential U(x) we may choose a w k = 
10 4 cm -1 . The number of atoms in the quantum soliton 
may fluctuate around 100, to be higher than the disin- 
tegration threshold, the latter being of the order of 1, 
according to the estimates given above. 

To conclude, we have revealed hitherto unknown prop- 
erties of quantum solitons traversing Gaussian potentials: 
the strong dependence of their transmission and disinte- 
gration upon the particle number. This property may al- 
low the observation of novel types of multipartite states. 

The support of the EC (the QUACS RTN and 
SCALA NOE), and ISF is acknowledged. I.E.M. also 
thanks the programs Russian Leading Scientific Schools 
(grant 1115.2003.2) and Universities of Russia (grant 
UR.01.01.287). 



[1] R.W. Boyd, Nonlinear Optics (Academic Press, NY, 
1992), §6.5; A. Newell and J. Moloney, Nonlinear Optics 
(Westview Press, Boulder, 2003), Chapter 6. 

[2] Y. Lai and H.A. Haus, Phys. Rev. A 40, 844 (1989); 40, 
854 (1989). 

[3] F.X. Karter and H.A. Haus, ibid., 48, 2361 (1993). 

[4] P.D. Drummond, K.V. Kheruntsyan, and H. He, J. Opt. 
B: Quant. Semiclass. Optics 1, 387 (1999). 

[5] R.K. Bullough and M. Wadati, J. Opt. B: Quant. Semi- 
class. Optics 6, S205 (2004). 

[6] R. Kanamoto, H. Saito, and M. Ueda, Phys. Rev. Lett. 
94, 090404 (2005). 

[7] S.A. Morgan, R.J. Ballagh, and K. Burnett, Phys. Rev. 
A 55, 4338 (1997); L. Khaykovich et al., Science 296, 
1290 (2002). 

[8] M. Machholm, C.J. Pethick, and H. Smith, Phys. 
Rev. A 67, 053613 (2003); M. Machholm, A. Nicolin, 
C.J. Pethick, and H. Smith, Phys. Rev. A 69, 043604 
(2004). 

[9] M. Greiner et al, Nature (London) 415, 39 (2002). 
[10] P.J.Y. Louis, E.A. Ostrovskaya, and Yu.S. Kivshar, Phys. 

Rev. A 71, 023612 (2005); 
[11] Th. Anker et al., Phys. Rev. Lett. 94, 020403 (2005). 
[12] H.A. Bethe, Z. Phys. 71, 205 (1931). 
[13] M. Girardeau, J. Math. Phys. 1, 516 (1960); E.H. Lieb 

and W. Liniger, Phys. Rev. 130, 1605 (1963); 

H.B. Thacker, Rev. Mod. Phys. 53, 253 (1981). 
[14] Z. Cheng and G. Kurizki, Phys. Rev. Lett. 75, 3430 

(1995); Phys. Rev. A 54, 3576 (1996). 
[15] LB. McGuire, J. Math. Phys. 5, 622 (1964). 
[16] L.D. Landau and E.M. Lifshitz, Quantum Mechanics 

(Oxford, Butterworth-Heinemann, 2000), §25. 
[17] Ph. Courteille et al., Phys. Rev. Lett. 81, 69 (1998); 

J.L. Roberts et al., ibid., 81, 5109 (1998); V. Vuletic, 



3 



